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Abstract 



rotating compact star in the gravitoelectromagnetic approximation to the 
geodesic equation . The star is assumed to be surrounded by an ideal, 
highly conducting plasma (taken as a magnetohydrodynamic fluid) with a 
' stationary, axially symmetric electromagnetic field. The general relativis- 

, tic Maxwell equations are solved to obtain the effects of the background 

spacetime on the electromagnetic field in the linearized Kerr spacetime. 
The equations of motion are then set up and solved numerically to in- 
corporate the gravitational as well as the electromagnetic effects. The 
Q ' analysis shows that in the slow rotation approximation the frame drag- 

^ I ging effects on the electromagnetic field are absent. However the particle 

C/3 . is directly effected by the rotating gravitational source such that close 

to the star the gravitational and electromagnetic field produce contrary 
effects on the particle trajectories. 



b ! 1 Introduction 



The study of accretion dynamics of charged particles in the vicinity of a com- 
pact gravitational source is of observational as well as theoretical interest; par- 
ticularly for the information it reveals about the nature and influence of the 
background spacetime on the various physical processes occurring in the star's 
vicinity. Whereas a compact star is formed largely of a degenerate plasma^, 
matter surrounding the star exists in a highly ionized, though a less dense, 
state^. In this region the charged particles are more free to move while not 
accreted to relativistic velocities, and the dynamics of these charged particles is 
governed not only by the gravitational but also by the electromagnetic field of 
the star. However the electromagnetic field is itself effected by the geometry of 
the background spacetime. Therefore it is of importance and of astrophysical 
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relevance that the two-fold effects; namely, the coupling of the gravitational 
and the electromagnetic field of the star, as well as the direct effects of the 
background spacetime on the motion of a charged particle, must be addressed 

together''*'^. 

Here we present an investigation of these effects on the dynamics of an ac- 
creted charged particle lying in the star's plasma atmosphere. To pose the 
problem quantitatively we consider a slowly rotating isolated compact star as 
the source of gravitational field which defines the geometry of the background 
spacetime. An axially symmetric, highly conducting plasma with a station- 
ary, axially symmetric electromagnetic field is assumed to constitute the star's 
atmosphere. We regard the electromagnetic field as perturbed by the back- 
ground spacetime, however the back-reaction of the electromagnetic field on the 
spacetime curvature is neglected and the standard Kerr metric (in its linearized 
form) is used for the fixed background spacetime. For the plasma surrounding 
the star the magnetohydrodynamic (MHD) approximation is assumed, as such 
the plasma is considered as an ideal highly conducting MHD fluid"''"^. The 
formulation of the electromagnetic field equations is based on a generalized def- 
inition of the electromagnetic field tensor for an ideal MHD fluid in a curved 
spacetimc^'^. The direct effects of the spacetime dragging on the particle dy- 
namics are discussed in the gravitoelectromagnetic (GEM) approximation to the 
geodesic equation^*^"^^. The GEM approximation, which assumes the charged 
particle motion to be slow, is generally valid in the case of stable compact stars 
(i.e. non-collapsing objects). Also for the case of a magnetized compact star 
with a dipole field the GEM approximation is applicable except at or very close 
to the poles where particle velocities can be relativistic due to intense magnetic 
field. Generally the accretion occurs well below the relativistic limit and so the 
GEM approximation can be used. Throughout we take the velocity four vector 
as defined by a comoving observer ZAMO (zero angular momentum observer). 
The ZAMO is a useful class of observers, relevant to astrophysics, circling the 
gravitational source with a given angular velocity at a fixed radial distance and 
polar angle^'^'^'*. 

The paper is organized as follows. In the first part of our analysis, comprising 
section H, we obtain exact solutions of the general relativistic Maxwell equa- 
tions using the generalized definition of the electromagnetic field tensor for an 
ideal, highly conducting MHD fiuid in linearized Kerr spacetime. The solution 
obtained shows that in the slow rotation approximation the effects of rotation 
on the electromagnetic field are absent. Then, in section III, assuming the star 
to be a sphere of homogeneous mass density, we study numerically the dynamics 
of the charged particle using the GEM approximation. In this approximation 
the direct effects of rotation (i.e. via frame-dragging) on particle trajectories 
are found to oppose that of the electromagnetic field, hence the particle motion 
is a result of balancing and counter-balancing of these effects. Lastly, in section 
IV, we give a discussion and a summary of the main conclusions of the paper. 
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2 The Electromagnetic Field 



2.1 General Relativistic Maxwell Equations 

The exterior spacetime to a slowly rotating compact gravitational source of 
mass M is described by the linearized Kerr metric. In the Boyer-Lindquist 
coordinates the metric can be written as the general line element: 



rfs2 = -e'^nr)(it'^_2uj{rysm^edtdip + r^sm'^edip^ + e-^^^^Ur^ + r^d0^, (1) 
where 

= (1 - (2) 

and 

^(r) = ^ = (3) 
dt 

is the angular velocity of a free falling frame brought into rotation by the frame 
dragging of the spacetime. Here and in what follows the Greek indices run as 
t, r, 6, and respectively. Also throughout we assume the gravitational units in 
which G = 1 = c. 

The general relativistic form of the Maxwell equations is 

-Fa/3,7 + -^/37,« + Pja,/} = 0, (4) 



(y^F"'^)^^, = 47rV^J", (5) 

where g represents the determinant of the metric tensor g^/s given by the Ein- 
stein field equations. Here Fa/} is the generalized electromagnetic field tensor 
for an ideal MHD fiuid given by a unique tensorial expression: 

Fa0 = UaEn - U^Ea + TJ^iSjsU'B^, (6) 

and J" is current four vector. In general the current four vector is the sum of 
two terms corresponding to a convection and to a conduction current: 

J« = eu" + au/sF/^" (7) 

where e is the proper charge density, a is the conductivity of the fiuid, and u" is 
the unit velocity 4-vector. The volume element 4-form rj^/s-ys and its dual r;"^'*"' 
are defined by 

Vai}j5 = V^ea^^s, r,"^''' = --^e"^''', (8) 

where ea^-ys is the Levi -Civita symbol, which is +1, —1, or for an even, odd, or 
non-permutation of afi-yS respectively. Also the four vectors Ea and Ba, denote 
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the electric and magnetic field components in the four dimensional spacetime. 
For a ZAMO Ur and u$ vanish, and using = — 1, the components of the 

four velocity vector are: 

= e-*('') (1,0,0, cj(r)), = e*W(-l, 0, 0, 0). (9) 

Here we employ the often made assumption of a perfectly conducting plasma 
surrounding the star (of e.g. [5], [14] for a discussion of this assumption with 
respect to the generalized Ohm's law). Since for plasma the condition of neu- 
trality also holds thus J" must vanish identically. The Maxwell equations (4) 
give: 

{^u'B-),r + {V^u'B'),e = 0, (10) 

{utE^).r = 0, (11) 

{utE^),e = 0, (12) 
{utEg - ^u'^B')^r - (utEr - ^u^B^)^e = 0, (13) 
whereas from expression (5) we obtain: 

{sf^u'E^)^r + {V^9u'E'),e = 0, (14) 

{utB^),r = 0, (15) 

{utB^),0 = 0, (16) 

{utBe - V^u'^E^)^r - {utBr - ^f^u^E'^)^e = 0. (17) 

2.2 Solution to the Mcixwell Equations 



The electromagnetic field outside the plasma smromiding the compact star is 
now determined by equation (10) to (17). To solve this system of equations let 
us assume the following separation ansatz for the magnetic field components B^ 
and B^: 

B^{r,9)=Ri{r)Q^j,{e), B'{r,9) = R%{r)Q%{0), (18) 
and similarly for the electric field components E'^ and E^: 

E^ir,e)=RUr)em, E\r,9) ^ R%{r)Q%{0). (19) 
Substituting from equation (18) into equation (10) and simplifying we obtain 

+ - — dTr — • ^2°) 

For separation we require that 0g = 8^ = 6b and obtain 
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Ob d6 
and 



1 '^^B +cot9 = ki, (21) 



^(^^"*^^)=fc„ (22) 



r^u*R% dr 



where ki is the separation constant. Solving for Qb the equation (21) we obtain 
Ai expki6{csc9). Notice that because of the exponential function the field at 
coincident points 6 = ■jt/2 and 7r/2 + 27r is not identical. We therefore require 
single valued solution ki = 0. The solutions to equations (21) and (22) are now 
given by: 



and 



^B = 4^, (23) 
sine* 



Rb = 4^, (24) 



where Ai and A2 are constants. Therefore we have 

B''{r,0)= , f. B\r,e) = R%{r)^, (25) 
where A = A1A2. And similarly from equation (13) and (18) we obtain 

sin^ 

and hence 



®E = 3z^ R'e = S^, (26) 



E'^{r.e)= ^ f . E\r,e) = R%{r)-?^ (27) 

where C = C1C2 and = 6^ = 0%. 

Now to determine the functions R^^ we substitute from (25) and (27) into 
equation (13). After some simplification we obtain in a separated form: 

d{utr^R%) _sme^^duj CcotO^ 

dr ~ Ci ^ dr r^sine'- ^ 

As before we have 



d{uyR%) 
dr 



= k2, (29) 



and 



A-^ sin e - — cos 6 - /C2C1 = 0, (30) 

dr r^ sm 9 
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where k2 is a constant. Since (30) holds for every value of 9, on comparing 
coefficients of sm.0, cos0, and unity we obtain C = = ^2 and A = 0. Here 
are now two possibilities for non-trivial solutions depending on either ^ or 
A^y^O (i.e. Ai = 0). 

For Ai equation (29) gives 

R% = A- (31) 

where A3 is a constant. Prom equations (25), (27) and (17) we obtain in a 
similar manner 



Hence 



and 



R'b = (32) 



B'{r,e) = Q, B'{r,e) = ^^ (33) 



E^{r,e) = Q, E\r,e)= ^f\ . (34) 

UtT'^ sm 6 



Also from equations (11), (12) and (15), (16) we obtain 



UtT^ Sin fj 



= , : 2. - (36) 



Utr'^ sin^ 9 

For Ai = (and similarly for Ci = 0), there is no poloidal component 
to the magnetic (electric) field and the only non- vanishing components to the 
electromagnetic field are given by equations (35) and (36). Clearly in both cases 
the electromagnetic field is independent of to{r). 

In Figures (1) and (2) we give plots for the magnetic (electric) field compo- 
nents (E^) and Bf (Ef) for A1A3 = 1 (= C1C3) and A4 = 1 {= C4). 

3 Trajectory of a Charged Particle 

According to the general theory of relativity the trajectory of a particle in the 
field of a massive object is determined by the geodesic equation: 

where F^^ is the Christoffel symbol. To linearize the geodesic equation we 
assume that the metric tensor can be expressed as: 

9oc0 = Vafi + hocf) (38) 
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where 77^^ = diag(— 1, — I, —1,1) is the Minkowski metric tensor and h^f} is a 
small perturbation to the spacetime metric such that hafj 1. We can take 
time t for the afRne parameter s. Further requiring v=\dr/dt |<C 1 for a 'slow' 
moving particle we obtain^^ 

J = (G + vxH) (39) 

where 



G = -V<^, H = Vx4a (40) 

and 

-IIIt"- 

p being the mass density and V is the volume of the gravitational source. This 
approximation, due to its formal analogy to the classical electromagnetic theory, 
is referred to as the gravitoelectromagnetic (GEM) approximation to the general 
theory of relativity (see reference [16] and [17] for different interpretations). 
Thus in expression (40) G is called the gravitoelectric (GE) force per unit mass 
whereas the term mv x H is called the gravitomagnetic (GM) force for a unit 
mass. In particular the GM force, being independent of the choice of a particular 
coordinate system used, is general relativistically significant. 

In the slow rotation approximation the deformation to the star due to the 
effects of rotation are small compared to the radially attractive GE force; hence 
we take the star to be a slowly rotating sphere of homogeneous mass density. 
Then using GEM analogy, we obtain: 

G = -Mr/r^, (42) 

and 

^ = ~MR'{n.v^-\^, (43) 

5 r' 6 r-' 

where M is the mass, R is the radius and Jl is the angular velocity of the 
star. Thus in the field of a slowly rotating compact star, the net force acting 
on a particle of mass m and charge q is the sum of the gravitational and the 
electromagnetic forces: 

=m(G + v X H)+g(E + v X B). (44) 

To specify the physical situation we study the motion of the charged particle 
in a Cartesian coordinate system (x,y, z) and take H to be along the positive 
2-axis. Furthermore, since now the effective component of the GM force lies in 
the Xy-plane, we assume that the XF-plane is the equatorial plane of the star 
and the particle's motion is confined to this plane only. 
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Then for the case of non-vanishing poloidal magnetic field we transform 
expressions (33)-(36) and (42)-(43) into the Cartesian coordinates, and obtain 
the equations of motion in the equatorial plane (i. e. for 6 = n/2): 



fx 



X"^ + 2/2 ^(a;2 + y2)2 _ 2M(x2 + 2/2)3/2 




+ 



( 



(a;2 + 2/2)3/2 ^^2_^^2_2Mv'a;2+2/2 



dy_ 
dt 



(45) 



di2 



+ y2 ^(2-2 _^ y2)2 _ 2M(a;2 + 2/2)3/2 



M 2/^^o 



( 




(46) 



(a;2 + 2/2)3/2 ^3.2 + y2 _ 2MVa;2 + 2/2 



where = gAiAs/m , Eq = qd/m and = 4Mi?2ri/5. Notice that here 
the case of vanishing poloidal magnetic field is obtained for Bo = whereas 
requiring Ci = does not alter the equations of motion (45) and (46). We 
numerically solve the equations of motion (45) and (46) for the initial conditions 
a;(0) = 10 = 2/(0) and dx/dt |t=o= = dy/dt \t=o, and plot (in Figure (3) to 
(6)) the solutions for various values of the parameters M, fi, Eo, and Bq. 

4 Discussion and Conclusions 

We have considered the dynamics of a charged particle in a highly, conducting 
plasma surrounding a slowly rotating compact gravitational source. The effects 
of background spacetime have been discussed, both directly on the particle's 
motion and indirectly via an axially symmetric electromagnetic field around the 
star. 

We found (equations (33) to (36)) that when slow rotation is assumed for 

the compact source the spacetime dragging does not modify or reduce the elec- 
tromagnetic field of the star. Hence the field remains, in this approximation, 
the same as for the Schwarzchild spacetime. However the motion of the particle 
is directly influenced via gravitational field of the star. 

Our numerical analysis of the effects of gravitational and electromagnetic 
field shows (Fig. (3) to (6)) that whereas the electromagnetic effects dominate 
the overall dynamics of the charged particle at sufficiently large distances, the 
gravitational field also becomes important closer to the surface of the star espe- 
cially if the star is rotating fast enough. Firstly we observe that the effects of 
the electric field is to lift a charge particle from close to the surface of the star 
(Fig. (3)) and make them move along the magnetic field lines. Here are now 
two possibilities: one, if the magnitude of the magnetic field is larger than the 
magnitude of the electric field the particle following a helical trajectory falls into 
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the star (case (a), (b), and (c) in Fig. (4)); two, if the magnetic field strength 
is weaker or is even comparable to the electric field strength, then the electrical 
effects become dominant and the particle escapes from falling into the star (case 
(d), and (c) in Fig. (4)). 

On the other hand the effects of gravitational field can be regarded as op- 
posing that of the electromagnetic field. In Figure (5) we notice that the GE 
force attracts the particle towards the star's surface against the electric field. 
Furthermore it is clear from Figure (6) that the gyroradius of the helical tra- 
jectory of the particle is increased due to the GM force; hence the GM force 
weakens the effects of the magnetic field. 

Summing up these observations we note that the accretion of a charged 
particle in vicinity of a compact star is mainly due to the GE and magnetic 
effects. These effects arc especially dominant at a sufficiently large distances 
from the star. However dose to the star the particle trajectory is also effected 
by the GM force as well as electric field which cause accerted charged particles 
to follow open lines of force. For a sufficiently fast rotating star and an electric 
field comparable to the magnetic field the above mechanism (GM -|- electric) 
may contribute to opposing the in-fall of the plasma surrounding the compact 
star. 
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Figure Captions 

1) Plot for the magnetic (electric) field component (E^) as a function of 
the radial distance r and polar angle 9 for A1A3 = 1(= C1C3). 

2) Plot for the magnetic (electric) field component B'^ {E"^) as a function of 
the radial distance r and polar angle 9 for = 1(= C4). 

3) Trajectory of the charged particle for electric force of varying strength 
with = 0.1, 1, 5, 10, 20, and M = 1, fi = 0.1, Bq = 10 in gravitational units. 

4) Trajectory of the charged particle for magnetic force of varying strength 
with Bq = 0, 10, 30, 50, 70, and M = 1, ^ — 0.1, _Eo = 1 in gravitational miits. 

5) Trajectory of the charged particle for GE force of varying strength with 
M = 1, 2, 3, 5, 10, and fi — 0.1, Eq — 1, Bq = 10 in gravitational units. 

6) Trajectory of the charged particle for GM force of varying strength with 
= 0.1, 0.2, 0.3, 0.5, 0.8, and M = 1, Eq = 1, Bq = 20 in gravitational units. 
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Figure 2 



Figure 3 




Figure 4 
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